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Проведена аналитическая оценка влияния центрального отражателя на примере плоского реактора в одногрупповом и диффу- 
зионно-возрастном приближениях. Со стороны внешнего отражателя вводится эффективная граница, что существенно упроща- 
ет математическую постановку критической задачи. Апробация решения проведена для одного из состояний исследовательско- 
го реактора ИРТ-Т. 


Введение 

Оценка критических размеров или критической 
концентрации (массы) горючего является основ- 
ной задачей теории ядерных реакторов. В обще- 
принятом изложении рассмотрение критической 
задачи ограничивалось реакторами классической 
геометрии (сфера, неограниченный цилиндр и 
пластина) без отражателя нейтронов и с его приме- 
нением с соблюдением центральной симметрии. 
Причем исторически сложившееся внешнее распо- 
ложение отражателя отвечало существующим на то 
время конструкциям ядерных реакторов. 

Влияние внутреннего отражателя рассматрива- 
лось только в одногрупповом приближении для 
предельных случаев его диффузионных характери- 
стик: для вакуума и черного тела. Такой подход по- 
зволял заметно упростить математическую поста- 
новку задачи (исключалось уравнение диффузии). 
Результаты решения критической задачи в одно- 
групповом приближении получили применение в 
расчетах отдельных типов реакторов (например, 
тяжеловодных реакторов) и в теории центрального 
поглощающего стержня. 

Более чем полувековой опыт эксплуатации 
ядерных реакторов выявил целесообразность при- 
менения центральных или внутренних отражате- 
лей. Например, в высокотемпературном реакторе с 
гелиевым теплоносителем ГТ-МГР отражатель 
представляет собой конструкцию, набранную из 
графитовых блоков различной конфигурации [1]. 
Он включает в себя не только обычные в примене- 
нии отражатели: верхний, нижний и боковой, но и 
центральный отражатель. 

Активная зона реактора ГТ-МГР состоит из 
призматических графитовых блоков и поэтому 
обладает меньшим термическим сопротивлением 
по сравнению с активной зоной из шаровых твэ- 
лов. Наличие центрального графитового отражате- 
ля повышает аккумулирующую способность реак- 
тора и тем самым снижает скорость разогрева ак- 
тивной зоны при ухудшении теплоотвода. Одно- 
временно такая кольцевая компоновка обеспечи- 
вает лучшее выравнивание энергораспределения 
по активной зоне и, как следствие, меньшие удель- 
ную мощность энерговыделения и температуру то- 
плива. 


В исследовательском реакторе ИРТ-Т примене- 
ние не только внешнего, но и внутреннего отража- 
телей обусловлено необходимостью увеличения 
средней плотности потока нейтронов [2, 3]. Этой 
же цели служит и выбор бериллия в качестве мате- 
риала отражателя. 

Распределение плотности потока тепловых 
нейтронов в общем случае проводится численны- 
ми методами в многогрупповом приближении и в 
многомерной постановке. Целью работы является 
получение аналитической оценки влияния цен- 
трального отражателя на примере реактора ИРТ-Т 
и с учетом принимаемых допущений носит скорее 
методический характер. 



Рис. 1 . Схема расположения элементов реактора ИРТ- Т 


Активная зона реактора ИРТ-Т в плане пред- 
ставляет собой прямоугольник со сторонами 6x8 
ячеек квадратного сечения. По периферии в один 
ряд располагаются бериллиевые блоки отражателя, 
в центре находится бериллиевая ловушка нейтро- 
нов (внутренний отражатель), включающая 4 ячей- 
ки, остальные ячейки занимают тепловыделяющие 
сборки, омываемые легководяным теплоносите- 
лем. В целом расположение элементов по вытяну- 
той стороне прямоугольника относительно оси 
симметрии имеет вид, представленный на рис. 1, 
где Ні/ 2 - границы соответствующих зон. 

Задачу о реакторе с внутренним и внешним 
отражателями будем решать в одногрупповом и 
диффузионно-возрастном приближениях, для про- 
стоты считая активную зону гомогенной и в форме 
неограниченной пластины. 
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Одногрупповое приближение 

Распределение потока нейтронов Ф, в активной 
зоне, внутреннем и внешнем отражателях в одно- 
групповом приближении принято описывать сле- 
дующими волновыми уравнениями: 

ДФ.+ЯІ 8 ®, =0, (1) 

АФ 2 -ХІ ф 2 = 0, (2) 

АФ 3 -Ж 3 2 Ф 3 = 0, (3) 

где Хъ Хъ Ъ — материальные параметры среды ак- 
тивной зоны, внутреннего и внешнего отражателей. 

Используем следующие граничные условия для 


решения уравнений (1-3): 
х = 0 ѴФ 2 (0) = 0, (4) 

х=Н І /2 Ф 2 (Я,/2) = Ф 1 (Я,/2), (5) 

х = Н 1 /2 0 2 ѴФ 2 (Я,/2) = 0 1 ѴФ 1 (Я і /2), (6) 

х = Н 2 / 2 Ф,(Я 2 /2) = Ф,(Я 2 /2), (7) 

х = Я 2 /2 ДѴФ 1 (Я 2 /2) = Я,ѴФ 3 (Я 1 /2), (8) 

х = Н 3 /2 Ф 3 (Я 3 / 2) = 0, (9) 

где Д - коэффициенты диффузии. 


Применение внешнего отражателя толщиной Г 
снижает утечки нейтронов и дает экономию актив- 
ной зоны равной эффективной добавке [4]: 

<5 эф =— агс{§ 1 ^-іЦХз Г > ■ 

Хі 1°ъХ 3 ] 

Упростим задачу ( 1-9), для чего введём эффек- 
тивное граничное условие 

х = Н 2 /2 + 5 эф = Я эф /2 Ф,(Я ф /2)=0. (10) 

Таким образом, задача на критическое состояние 
будет определяться системой ур. (1, 2, 4-6) и (10). 

Решения волновых уравнений (1) и (2) извест- 
ны и имеют вид: 

Ф,(х) = С, 8 ІПІХіХ) + С 2 со8(^х), 

Ф 2 (х) = С 3 к1і(х 2 х) + С 4 сй(/ 2 х). 

Удовлетворяя эти решения граничным усло- 
виям, получим систему алгебраических уравнений. 
После снижения порядка системы до двух и затем, 
освободившись от констант, получим критическое 
уравнение для реактора в форме неограниченной 
пластины с центральным отражателем: 

О 2 *2 1| і( Х 2 я і /2) = -Я, ^сІ8(^,(Яф /2 - Я, /2)). ( 1 1 ) 

Из этого уравнения при заданной толщине вну- 
треннего отражателя и известных материальных 
характеристиках можно найти критическое значе- 
ние размера активной зоны. 

Исходная система алгебраических уравнений 
позволяет найти выражения пространственного 
распределения плотности потока нейтронов толь- 
ко с точностью до какой-либо константы С,. Для 


нахождения законов распределения нейтронов в 
активной зоне и отражателе в явном виде зададим 
дополнительное условие: 


х = 0 Ф 2 (0) = Ф* = длд (12) 


и в результате получим законы распределения [6]: 


Ф 1 (х)=Ф*сЬ(х 2 Я 1 /2) 


5ш(Хі(Яз/2-х)) 

8іп {х,{Н г І2-Н 1 /2)У 


(13) 


Ф 2 (х)=Ф*сЬ(х,х). (14) 

Если принять, что Ф*=1, тогда решения (13, 14) 
будут иметь вид, представленный на рис. 2. 



О 10 20 30 х, см 

Рис. 2. Распределение потока нейтронов по активной зоне и 
центральному отражателю 


Двухгрупповое приближение 

Более высокое диффузионно-возрастное при- 
ближение учитывает замедляющиеся быстрые и 
диффундирующие тепловые нейтроны. Также вос- 
пользуемся эффективным граничным условием. 



В соответствии со схемой двухгрупповых про- 
цессов (рис. 3) критическая задача в осесимме- 
тричной декартовой постановке будет включать 
следующие уравнения баланса нейтронов для: 

• активной зоны Н,/2<х<Я,/2: 
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ОАФ І -Ъ к Ф [+ к а ? аТ Ф 2 / ( р= 0, (15) 

ДЛФ 2 -Е оД .Ф 2 4-фІ^Ф, = 0; (16) 

• внутреннего отражателя 0<х<#,/2: 

А,ЛФ„.-^,Ф 1( = О, (17) 

А,АФт,-2„, г ,.Ф 2 ,.+2 й ,Ф, =0. (18) 

Граничные условия представим в виде: 
х = 0 ЛФ 1( , = 0, АФ 2г = 0, (19) 

х = Н 1 / 2 Ф,ДЯ,/2) = Ф,(Я,/2), 

А,ѴФ 1г (Я 1 /2) = ДѴФ 1 (Я 1 /2), (20) 

х = Н 1 /2 Ф 2г (Я 1 /2) = Ф 2 (Я, /2), 


А, (А /2) = А ѴФ 2 (А /2), (21) 

х = Н 2 / 2 Ф,(Я,/2) = 0, Ф 2 (Я 3 /2) =0, (22) 

где 2Д Г , - макроскопические сечения поглоще- 
ния тепловых и увода быстрых нейтронов, см -1 ; А, 
(р - коэффициенты размножения бесконечной сре- 
ды и вероятности избежать резонансного захвата. 

Законы распределения потоков нейтронов в актив- 
ной зоне и отражателе 

Поскольку уравнения (15) и (16) симметричны 
относительно потоков, то это дает основание счи- 
тать, что потоки оцениваются волновыми уравне- 
ниями с одним и тем же волновым числом [4, 5(: 

АФ, + а 2 Ф, = 0, ДФ,+сгФ 2 =0, (23) 

откуда выразим лапласианы в виде 

АФ, = -сГФ,, АФ, = -а 2 Ф 2 
и подставим их в (15) и (16). В результате получим: 
(а А>, +5Д)Ф, = АА/Ф 2 /ф> 

(а 2 А+А,г)Ф 2 = <РАФ,- 

После деления одного уравнения на другое и 
введения квадрата длины диффузии тепловых ней- 
тронов І 2 = Д/Е а г и их возраста т=Д/2Д, после не- 
сложных преобразований получим: 

к., = (1 + а А 2 )-(1 + а 2 г). 

Это квадратное уравнение относительно а 1 
имеет два корня: а, 2 =/Г и а 2 2 =-ѵ 2 , где А и ѵ 2 - чи- 
сленные константы, аналогичные материальному 
параметру [4] : 



Таким образом, ур. (23) для быстрых нейтронов 
распадается на уравнения: 

ДФ, + ц 2 Ф, = 0, АФ, -ѵ 2 Ф, = 0 

и его решение есть линейная комбинация решений 
последних уравнений: 

Ф, (х) = А,8іп(цх) + А 2 со5(/лх) + 

+А 3 зЬ(ѵ х) + А 4 сЬ(ѵ х). (25) 

Аналогично получим общий вид решения для 
потока тепловых нейтронов: 

Ф,(х) = 5,5іп(цх) + 5 2 со8(цх) + 

+5,8іі(ѵх) + й 4 с1і(ѵх), (26) 

где Д=5Д„ а 5, - соответствующие коэффициенты 
связи. 

В уравнении диффузии для быстрых нейтронов в 
отражателе ( 17) отсутствует член генерации быстрых 
нейтронов, поскольку среда не мультиплицирую- 
щая, поэтому система уравнений (17, 18) перестает 
быть симметричной и распадается на два последова- 
тельно решаемых уравнения. Однородное уравне- 
ние для быстрых нейтронов, запишем в виде: 

АФ -В 2 Ф =0 

где р (ДЕдДѴ Общее решение этого уравнения из- 
вестно, и его можно записать в виде: 

ф і,А) = Д -8ЦР 1г х) + Р 2 •сіі(Д.х). 

Поскольку должно выполняться условие симме- 
трии (19), отсюда следует і 7 ,=0, и решение примет вид: 

ф ъ А) = А • сЬ( р 1г х) = Р ■ сЬ(Д. А- ( 27 ) 

Поток тепловых нейтронов в отражателе опре- 
деляется ур. (18). Найдем вначале решение одно- 
родной части этого уравнения: 

А г АС-АД:=о 

или с вводом параметра Р\=У. аТг / О г =\ / П г 

А ф І-РІ ф І = 0. 

Аналогично, в силу условия (19), его решение 
будет определяться зависимостью: 

Ф 2 °г А) = А • сЬ( Р 2г х) = С ■ с1і( р ъ . х), 

а решение неоднородного уравнения (18) можно 
записать в виде: 

Ф 2г А) = С • с1і(/3 2 ,.х) + 8 Г Р ■ сЫР ь . х). (28) 

Для удобства общие решения уравнений 
(15-18) перепишем в виде: 

Ф і (х) = А 1 Х і +А 2 Х 2 +А 3 У 1 +А 4 У 2 , 

Ф-) А) = А, ДХ, + 8 п А 2 Х 2 + 8 3 А 3 У 1 + 8 4 А 4 У, , 

> (/У) 

Ф 1г (х) = РХ 1г , 

Ф 2 Хх) = ОХ 2г + 8РХ ]г , 

где Д, Р, С - константы, определяемые из гранич- 
ных условий; функции - Х,=$іп(,их), Х 2 =со $(/их), 
Г|=8Іг(ѵх), Г 2 =с1і(ѵх), X, =с1і(Д„х), Х 2 =с ЦР 2 рсУ, А А 
А А Д - коэффициенты связи. 
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Поиск коэффициентов связи 

Процедура определения коэффициентов связи 
подробно изложена в [4, 5], покажем её на примере 
нахождения значения Д. Подставим неоднородную 
часть решений Ф ]г (х)=РХ ]г и Ф }г (х)=Х г ГХ ]г в ур. (18): 

ДД™ 1г -і л ^,н 1г +Д.ЛЛ, = о. 

С учетом уравнения А Х ІГ =Х ІГ /т получим 

В гА рх ь ./т - Д Тг Д РХ ь ., +Х Кг ГХ и . =0, 


Детерминант системы (34) при заданной толщи- 
не внутреннего отражателя и известных материаль- 
ных характеристиках (условие критичности) позво- 
ляет найти эффективное значение критического 
размера активной зоны. Экономия отражателя 
определяется по известной формуле [4] и позволяет 
найти истинное значение критического размера. 

Определение констант в уравнениях 
распределения потоков нейтронов 


откуда 


А =- 




X ... — О, / т 

а,Тг 2 г г 


(30) 


: Т (т -Ь 2 )' 

а,Тг\*г г/ 

Аналогично определяют другие коэффициенты 
связи: 

_ <Р^я 


А,=5 2 = 


Х„ г +(х ] В 2 Е в7 (1+а,Х-) 


(31) 


Д=^ = 




Е г (1-а 2 2 ! 2 )' 


(32) 


Выражения (30-32) для коэффициентов связи 
остались прежними [4, 5], что объясняется иден- 
тичностью исходных уравнений. 


Условие критичности 


Для определения констант в общих решениях 
уравнений (29) удовлетворим их граничным усло- 
виям (20—22): 

А 1 [Х 1 ] + А 2 [Х 2 ] + А 3 [У 1 ] + А 4 [Ѵ 2 ]=Г[Х іг і 

дд[хд+дэ 2 [х 2 ]+дэ 3 [д]+дл 4 [г 2 ] = 

= ^Г[Х и .] + С[Х 2г ], 

ДЛ , [' ѴХ, ] + Д А 2 [ ѴХ 2 ] + ДЛ з Г ѴД ] + 


+ДЭ 4 [ѴГ 2 ] = ДЛѴД,.], 

А, Д А , [ѴД ] + ДДЛ 2 [ѴХ 2 ] + Д Д Э , [ ѴД ] + 
+ДДЭДѴД] = ДДД[ѴД.] + Д,.С[ѴЛ ], 
4 1 Д 1 ) + 4 2 (1 2 ) + 4 3 (Д + Л(Г 2 ) = 0, 


(33) 


ДД(ХД + ДЛ 2 (Х 2 ) + ДД(Д) + ДЛ(Д) = 0.] 


Напомним, что граничные условия (19) мы уже 
использовали, что позволило получить уточненные 
решения (27) и (28). 

Квадратные скобки использованы здесь для 
обозначения функций и их производных на грани- 
це активной зоны и отражателя, круглые - на эф- 
фективной границе зоны. 

Система (33) имеет нетривиальное решение 
только в случае, если ее определитель равен нулю: 
д = 

га га пл щ га ] о 

д,га д,га зді з, ру .ѵ,га: -гад 

Ц[ѴТГ,] Ц[ѴА] Ц[Ѵіу Ц[ѴУД -ЦДѴХ,,] о 

_ Д,А[ѴД’,] дддѵтг,] д,А[ѴТ] ДАІѵгу -ААігаі -АДѴА,] 

(тг,) (А) А) (К) о о 

5 , га д,(А) ЗД) ад) о о 

= о. 


(34) 


Чтобы вычислить распределение потоков в ак- 
тивной зоне и отражателе необходимо знать крити- 
ческий размер активной зоны и постоянные А ь А ъ 
А ь А 4 , Си Г, ур. (29). Поскольку число неизвестных 
превышает число исходных уравнений, то законы 
распределения можно найти лишь в неявном виде 
с точностью до определенной константы. Для уста- 
новления законов распределения нейтронов по ак- 
тивной зоне и внутреннему отражателю в явном 
виде воспользуемся дополнительным условием: 

Х = 0 Ф 1г (0) = ФГ, 


откуда следует 

ФГ = С- сЬ(Д.0)= С. 
Другие константы примут вид: 


А, = 


ФЛДД-Л[Х 2 ]-ДД]-Э 4 [У 2 ] 


[Д] 


-ф п 


А 0 =- 


Лдд- ' 

- д 

'т- 

{-р[хх 1г ]) 




А 4 [ѴГ 2 ](г а[У 2 ]) 


А _ -фГ8 1 ~аа 


[ѴХ 2 ]-а[Х 2 ] 

ФГ (8А-8&) 


А,= 




<5 , = Р[Х Ь .]- ([ѴХ 2 ]-а[Х 2 ])~ (а[Х ѵ ] - 
-Р[ѴДД)-((Х 2 )-Д[Х 2 ]), 
5 2 =(«)-ДК])-([ѴХ 2 ]-а[Х 2 ])- 
-([ѴД] -а[Д])-((Х 2 )-Д[Х 2 ]), 
5 3 =((Д)-ДК])Ч[ѴХ 2 ]-а[Х 2 ])- 
-([ѴГ 2 ]-а[Г 2 ])-((Х 2 )-Д[Х 2 ]), 
Д=ДД[Д,.]([ѴХ 2 ]-а[Х 2 ])- 
-(а[ДД- р[ѴД,,])- (Д(Х 2 )- ДД[Х 2 ]), 

Д = (Д (Д ) - Д р [ Д ]) • ([ѴХ 2 ] -а[Х 2 ]) - 
-([V Д ] - « [ Д ]) • (Д (АС, ) - ^ >3[ аг 2 ]), 

Д = (Д (У 2 ) - Д Р [Г 2 ]) • ([ѴХ 2 ] - а[Х 2 ]) - 
-([ѴГ 2 ]-а[Г 2 ])-(Д(Х 2 )-ДД[Х 2 ]), 
Р=д г /А. а=[ѴХ 1 ]/[Х 1 ], Д=(Х 1 )/[Х 1 ]. 


Обсуждение результатов 

Проверка решения сделана на примере реакто- 
ра ИРТ-Т, активная зона которого набрана из вось- 
митрубных тепловыделяющих сборок. На конец 
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кампании на мощности 12 МВт реактор имеет сле- 
дующие характеристики активной зоны: 
г=45,831 см 2 , Д=1,632 см, Д=0,292 см, 
2 аГ =0, 10352 см' 1 , <р=0,973, ^=1,3463 и отражателя: 
т= 92,037 см 2 , 7) 1г =0,538 см, А г =0,373 см, 
2 аГ =0, 00873 см -1 , І 2 =42,745 см 2 . При построении 
графиков распределения нейтронного потока удоб- 
но принять Фі™ п =1, и тогда для нашего случая коэф- 
фициенты будут равны: 4=0,472388, М г =2,528782, 
4=85,387631, А=— 85,387631, Р=\, №=-0,57191. 
Для определения среднего потока тепловых ней- 
тронов по активной зоне необходимо найти эффек- 
тивную добавку отражателя. Результаты решения в 
среде МаДіСАЕ) представлены на рис. 4. 

Эффективность отражателей принято оцени- 
вать не только по экономии активной зоны, но и по 
изменению отношения среднего потока тепловых 
нейтронов к его максимальному значению. Для ре- 
актора без отражателей оно известно [4, 5] и равно 
0,637. Если применить внешний бериллиевый 
отражатель, то это отношение составит 0,939. 

В случае применения центрального отражателя 
(рис. 2) значение этого отношения в одногруппо- 
вом приближении по данным [6] получилось рав- 
ным 0,82. Этот результат не является неожиданным 
и может быть объяснен в первую очередь тем, что 
одногрупповая модель не учитывает вклад в нако- 
пление тепловых нейтронов в отражателе за счет 
замедления быстрых нейтронов. 
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Решения в двухгрупповом приближении 
(рис. 4) наглядно показывают роль центрального 
отражателя. Эффективность нейтронной ловушки 
позволяет поднять отношение среднего потока те- 
пловых нейтронов к его максимальному значению 
до 0,966. 



Следует отметить, что при переходе к эффек- 
тивной границе деформируется распределение 
нейтронного поля вблизи границы раздела зона - 
внешний отражатель. Более детальный учет ра- 
спределения потребует увеличения порядка систе- 
мы и заметного усложнения выкладок. 
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